Weak Massive Gravity 



o 

PL, 
OO 



Ph. 
43 



> 

O 

m 



x 

S3 



D. Comelli a , F. Nesti b and L. Pilo 6 ' c 
a INFN - Sezione di Ferrara, 1-35131 Ferrara, Italy 
b Dipartimento di Fisica, Universita di L'Aquila, 1-67010 L'Aquila 
C INFN, Laboratori Nazionali del Gran Sasso, 1-67010 Assergi, Italy 
(Dated: February 20, 2013) 

We find a new class of theories of massive gravity with five propagating degrees of freedom where only 
rotations are preserved. Our results are based on a non-perturbative and background-independent 
Hamiltonian analysis. In these theories the weak field approximation is well behaved and the static 
gravitational potential is typically screened a la Yukawa at large distances, while at short distances 
no vDVZ discontinuity is found and there is no need to rely on nonlinear effects to pass the solar 
system tests. Fhe effective field theory analysis shows that the ultraviolet cutoff is (mMpi) 1 ^ 2 ~ 
l//xm, the highest possible. Fhus, these theories can be studied in weak-field regime at all the 
phenomenologically interesting scales, and are candidates for a calculable large-distance modified 
gravity. 



Introduction. The long standing quest for a healthy 
modification of gravity at large distances has recently 
been the target of renewed interest. Since the work of 
Fierz and Pauli []| (FP) it was realized that in generic 
Lorcntz-invariant (LI) massive deformations of gravity 
6 DoF are present already at linearized level around 
Minkowski space, the sixth mode being a ghost. While a 
tuning can get rid of this mode at linear level, it reappears 
at the nonlinear level or around non-flat backgrounds ■ 
It is remarkable that only recently a nonlinear comple- 
tion was found which propagates five DoF also in gen- 
eral. This was proposed first up to the fourth order in 
perturbation theory 0] and then to the full nonlinear 
level in 0, @ • At linearized level this theory reduces to 
the Lorentz-invariant FP theory, and shares with it the 
failure to reproduce the correct light bending, at odds 
with General Relativity (GR) and observations. The is- 
sue does not disappear in the limit of vanishing graviton 
mass, the so called vDVZ discontinuity [(|. A way out 
was proposed by Vainshtein in 0, by arguing that non- 
linear effects restore the correct GR behavior at short 
distances from a source. While this mechanism has been 
verified in a number of models § , the theory has to rely 
on strong nonlinearities even at the macroscopic solar 
system scales, where the gravitational potential is small, 
as discussed for generic LI theories in Q. 

A different perspective is provided by more general 
theories of massive gravity where Lorentz invariance in 
the gravitational sector is not imposed, though the weak 
equivalence principle still holds. At the quadratic level a 
number of such theories that are free from the vDVZ dis- 
continuity is on the market EMI- 

Moreover, in some 

of them the cutoff beyond which the theory cannot be 
trusted or where nonlinearities outbreak, has been ar- 
gued to be the highest possible [Icj. Clearly, if there 
exists nonlinear completions of such theories, they rep- 
resent serious candidates for a consistent and calculable 
large distance modifications of gravity. 

In this Letter we uncover a full class of nonlinear the- 
ories implementing these ideas. The analysis is based on 
a recent work [l5| stating the general conditions for the 



propagation of five DoF in massive gravity, whose main 
result is reviewed in the next section. In the following 
section we provide the general solution to these condi- 
tions. In the remaining part we describe an explicit ex- 
ample providing massive gravity theories which are free 
from the vDVZ discontinuity at short distances, lead to 
screened (Yukawa) gravitational potential at large dis- 
tances, propagate five massive DoF both at linear and 
nonlinear level, and are calculable in the weak-field limit 
at phenomenologically interesting distances. 

Massive gravities with 5 DoF. The generic La- 
grangian for massive gravity is obtained by adding a non- 
derivative potential V for the metric g^, 
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J d 3 x^g- 
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V(g) + Cr, 



(1) 



and its propagating DoF were analyzed in [la] using the 
Hamiltonian formalism, with the standard ADM decom- 
position [lj| in terms of lapse TV, shifts TV 1 and spatial 
metric jij 

'-N 2 + N i N^ ij y i:j N^ 
7y -TW 7ij - 



9^ 



(2) 



The potential V is regarded as a function of these vari- 
ables, and it is useful to define V = m 2 TV det^) 1 / 2 ^. 
We will also denote N A = (TV, TV*) with A = 0,1,2,3 
(note TV is not a Lorentz 4-vector), and with Fab... the 
derivative of any F with respect to N A , N B , etc. 

As it stands, the presence of V in the action violates 
diffeomorphism invariance, and in principle six modes 
can propagate, with a (sixth) mode being typically a 
ghost. The conditions on the potential V(TV^, j* 3 ) to 
have 5 propagating DoF and unbroken rotations were 
found in [la ], and they turned out to be quite simple: 
• First, the four secondary constraints have to determine 
only three out of four auxiliary fields [Til ] , leaving a 
new constraint. For this, the hessian of the potential 
with respect to lapse and shifts must have rank 3 



det Vab — , (rank Vab = 3) . 
So, Vab has only one null eigenvector \ ■ Vab X 1 



(3) 
= 0. 
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last auxiliary field [lg], to complete the new constraint 
in phase space. This amounts to the following additional 
condition 



which can be now readily solved by stating that U is a 
scalar function of a particular combination of £ l and y l ° : 



U=U{X 13 ) with X i 



IT 



(8) 



= 0. 



(4) 



Here we defined V = 7~ 1 / 2 V and £ A = x A /x°- We are 
free to normalize to x = 1. 

It can be shown that the same eqs. ©, (01 hold also for 
explicitly time dependent potentials [18j . One can verify 
that the Lorentz-invariant massive gravity theory of [3, [B[ 
nontrivially satisfies both (J3|) and @. However, we show 
here that a wider class of solutions with phenomenologi- 
cally interesting features exist when the potential is not 
Lorentz invariant but only rotational invariant. To this 
aim, let us the discuss the general solution to ((3]) and (j4} . 

General solution. We start from eq. ([3]), by recalling 
that it is a 'homogeneous Monge- Ampere' equation, a 
well studied problem in mathematical physics. In the 
very clear work [l9j . the general implicit solution was 
found by changing variables from N A to the eigenvector 
\ A ■ In these variables, one finds that the potential must 
be an homogeneous function of \ A of degree zero. Then 
the key role is played by a generic function U(Cjij), 
defined as U — \ A Va and linked to Vo by a Legcndre 
transformation 



Vi = Ut 



(5) 



Note that V A = dV/dN A while K t = 8U/dg. Any func- 
tion U will lead to a solution of the Monge- Ampere equa- 
tion (J3J), provided dp can be integrated to find V. Inte- 
grability is ensured [19( by an important implicit relation 
between N l , N and £ l p which includes a new generic 
function £(f , y^): 



N l - N C = U" £ 



3 ' 



(6) 



where W 3 = (Uij)^ 1 (inverse of the hessian matrix) and 
d = d^iC. For any C, one should formally invert (JB]) to 
determine C(N A ,y ij ), and then integrate © in N°, N { 
to find the potential. As a result, the general solution 
of §5§ is expressed in terms of two generic functions U 
and C of £ 4 and 7^ . Note that up to now the variables 
jij inside U and C are spectators. 

Equation (|4]) brings jij into the game. When written 
in terms of U, it reads 



dU_ 



2^ 



dU 

dyV 



(7) 



1 We note en passant that this implicit change of variables be- 
tween TV 1 and £ l trivializes to a good extent the complicated 
Hamiltonian analysis of Il5ll . and contain the shift transforma- 
tions of Q • This points to a possible important physical role of 
the £* variables, which in ADM language are linked to a 'spatial 
over temporal' speed of the foliation. 



The last step is to find the potential by integrating the 
relation ([5]). These equations admit, as integration con- 
stant, a generic function C(7 U ) to be added to V. 

Summarizing, we have found a whole new class of so- 
lutions propagating five DoF, which are parametrized by 
three generic functions U{X % i), £(£*,7 y ) and C(7 y ). In 
view of unbroken rotations, IA has to be a scalar function 
of X lk 8kj, i.e. of its three invariant traces, and similarly 
for C = C(Y k S kj ). 

In order to integrate ([5]) and find a potential V ex- 
plicitly, some concrete C has to be chosen first, to solve 
the implicit relation © between and N A . We remark 
that in general the resulting potential, function of N, N l 
and "fij, will break Lorentz symmetry. The conditions 
for Lorentz invariance, leading to the potentials of [3, [f|, 
consist in complicated additional constraints that will not 
be discussed here. Instead, interesting new solutions can 
be found by many choices of C. The general analysis will 
be reported in a separate publication [l8| — here, it is 
enough to consider the case C = 0, which already leads 
to interesting features and phenomenology. 

An example of explicit theories. For C = 0, cq. (|6j 

leads to the simple relation = N' l /N. This gives 
XV = (Y 3 - N' l N J ' /N 2 ) = gv, which remarkably co- 
incides with the spatial part of the inverse metric. It is 
pleasing that also equations ([5]) can be integrated and 
one finds the following large family of explicit potentials 



V = U{g ik 5 kj )+N- 1 C(y ik 5 kj ). 



(9) 



It is interesting that in this example, a closed form for 
the Hamiltonian of the system upon using the secon dary 
constraints can be found. Using the expressions in [la ], 
we find that U cancels out and (suppressing boundary 
terms) 



H 



Ml m 2 



d 3 Xy/jC. 



(10) 



As a result, by choosing C a positive function, the Hamil- 
tonian is positive definite. 

A basic test of the above theories consists in analyzing 
the weak field expansion around flat space, to see whether 
it is healthy and the phenomenological implications are 
consistent with the weak field tests of gravity. 

The conditions on the potential ©, to admit Flat 
Minkowski background rj^ v are 



U \g =r] 



0, 



0, 



(11) 



\g=v 



where \ g=n = {N ->• 1, A' 0,7y -> <%} and d^C 
C'Sij, dyijU = U'Sij on Minkowski. 
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At quadratic level, the expansion of the potential leads 
to mass terms for the metric fluctuations = Qu,v~f)u,v> 
which can be parametrized as 

j r2 

pl I 2 7.2 i o 2 t.2 2 i 2 , 2 i2 o 2 j j\ 

— — yn h 0Q + 2m 1 h 0l — m 2 + m 3 h — 2 m 4 h o n j 

(12) 

with h = hu and repeated spatial indices are summed. 
The quadratic expansion reflects the breaking of Lorentz 
invariance in the full potential and gives the most general 
rotationally invariant mass terms for the gravitono They 
were studied in [13, [ll| , where it was shown that absence 
of ghosts and stability of the linearized theory can be 
achieved for m = 0, together with the conditions m\ > 
and ml > m\ > 0. It is thus crucial to verify whether 
the potential © can satisfy these criteria. Expanding 
the potential ©, we get the mass spectrum 

ml =2m\= m'\ g=n , ml = 0„. u \ g=j] (13) 

Oiiiuva* = d i^H ( d imi V - \ 7<i3i V ) • 

We thus find a strong correlation between m\ and 7714, 
and their ratio m\jm\ = 2 is nicely consistent with the 
stability conditions mentioned above0 Concerning the 
conditions m\ 4 > 0, the class of potentials is large and 
they can generically be satisfied. 

As expected from the non-perturbative canonical anal- 
ysis, conditions ([3]) and (fj| to have 5 DoFs have led to 
m\ = 0, which implies the absence of the ghost sixth 
mode at the linear level. Needless to say, after the anal- 
ysis of [IH , in these theories five DoF propagate both at 
linear and at nonlinear level. Notice the important role 
of the function C, i.e. the integration constant of (|5|): if 
C = then m\ = and only 3 DoF propagate at the 
quadratic level, i.e. the theory is strongly coupled at all 
scales. 

Still at quadratic level, stability of LB theories was 
studied also in FRW curved spacetime in For in- 

stance in deSitter space with curvature H 2 , five stable 
modes propagate provided H 2 < m\/{Z{m\ — to§) + to 2 .). 
This means that the graviton mass scale has to be larger 
than the Universe curvature, or in other words that it 
can take physically interesting values. 

Stiickelberg Point of View. In any theory with a 
broken gauge symmetry, the symmetry can be recovered 
by introducing extra Stiickelberg fields. In our case, in 
the unitary gauge, the potential (j9)) depends on the ro- 
tationally invariant scalars made out of the 3D matrices 
X* = g lk Skj, = -f lk Skj, as well as on the 3D scalar 
N. The breaking of diffeomorphisms can be understood 
as due to the explicit presence of N and of the external 



2 The Lorentz- Invariant FP case is recovered for mo = 0, rrn = m. 

3 Note, this relation is modified in the general case C ^ 0, and the 
stability conditions turn into conditions on C llgll . 



'frozen' metric Sij in the potential, and can be restored 
by using four Stuckclbcrg fields $ and for instance 
by promoting X, T and N to manifestly gauge invariant 
quantities (see (To, [llj for a thorough analysis). 

In the standard expansion around Minkowski space, 
<1> = t + <j>°, = x % + <j) 1 , the metric fluctuations in 
the quadratic Lagrangian (TT2l) arc replaced by the gauge 
invariant combinations = h^ v — 2<9( /i </>„) . Because 
TOq vanishes, <j)° does not propagate. On the contrary 
the three spatial Stiickelberg fields do. One can also 
decompose them in transverse and longitudinal parts, 
<t>i = At + diip, which makes clear the difference with 
the known Lorentz-invariant case. Indeed, while in the 
Lorentz-invariant case the longitudinal field ip is made 
canonical through rescaling by l/m 2 M p , which leads to 
low strong coupling scale discontinuity and Vainshtein 
radius, here all three fields are made canonical by the 
minimal normalization ~ l/mM p . In particular for tp 
this fact is due to the difference m\ — m\ being nonzero 
(see (fl4|) below) and is thus a strict consequence of 
Lorentz-breaking in the gravitational sector. 

Together with the standard canonization of the gravi- 
ton this normalization leads to three important conse- 
quences in the full theory: i) the cutoff of the theory is 
A2 = [m Mpi) 1 / 2 [l(|; ii) via mixing with the metric, the 
coupling of the Stiickelberg fields to matter sources van- 
ishes as to, and thus no discontinuity is present; iii) again 
in the presence of sources, there is no scale of nonlinear- 
ity (Vainshtein) larger than the standard Schwarzschild 
radius R s (provided A 2 > l/R s )- 

It is also instructive to show the quadratic Lagrangian 
of canonical propagating fields in the 'decoupling limit', 
defined as to — > 0, Mpi — > 00 keeping A 2 = \J to Mpi and 
T/Mpi finite. One has symbolically 

c eff = h c a hc + (^) 2 - ^iK -^i) (v )2 + 

2m\{m\-ml) 

+ (AD 2 -^|(V,A C ) 2 + ^ (14) 
m 1 Mpi 

which has two transverse tensor DoF in h c , two in the 
transverse vector A\ and one in the scalar ip c . All these 
modes are massive, with different masses (they is not 
visible in the decoupling limit, but see [3). Note that 
the scalar kinetic term in ([T4| is singular in the LI limit, 
where all -> m. Lorentz-breaking is present in the 
dispersion relations of vector and scalar, which can prop- 
agate at speeds different from the one of the light. The 
presence of these three extra polarizations, in addition to 
the interesting massive nature of transverse gravitational 
waves (see e.g. [22j), could thus lead to striking signals 
already in the forthcoming gravitational waves searches. 

Phenomenology. Let us focus on the gravitatiaonal 
potential generated by a static source, with a non-zero 
energy- momentum-tensor component too- 

At linearized level, the gauge invariant gravitational 
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potential $ = hoo/2 is (we set Mp = 1 here) 



A 



4> 



-3™3 
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2 A 2 
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M| 
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2 ' 



'llll 



(15) 



where the squared masses -M 2 2 ~ to 2 an d the dimen- 
sionless coefficients A± t 2 in the formal decomposition of 
the second line depend on the mass parameters m 2 , and 
A± + A2 = 1 [13] . The squared masses M 2 2 may turn 
out to be real or complex. As a result, the radial profile 
of the gravitational potential from a point source can be 
Yukawa, oscillatory, damped (or even plain exponential 
or Newton, for some special choices of the to 2 .) Once 
the squared masses M 2 and M| are real positive, the 
linearized potential is a sum of two Yukawa terms: 



M 



$ = — (A 



1 e 



-Mir 



.4 



2 e 



-M 2 r 



(16) 



At small r both gauge invariant potentials [T(| reduce to 
the linearized GR result, i.e. there is no discontinuity. 

In the full nonlinear theories that we propose, we can 
go one step further. In fact, the crucial difference with 
the Lorentz invariant case is that here we can use pertur- 
bation theory. This is actually required to confront with 
standard post-Newtonian (PN) tests. Using second order 
perturbation theory, we find that the vacuum spherically 
symmetric solution can be written, in appropriate New- 
tonian coordinates (r, p), in the form 



ds 1 = - 



2M 2 2 / M M 2 \ 

1 h m p ci h c 2 — 5- mimp) 

P \ P P J 



dr 2 



P 2 dn 2 



2M 2 J M M 2 
1 h to p C3 h C4— 5- 

p \ p p z 



dp 2 (17) 



where for convenience the metric has also been expanded 
to leading order in mp -C 1. The constants Cj depend on 
the chosen potentials U and C. This result shows that 
the perturbation expansion near a source has no Vains- 
htein cutoff, is valid as in GR down to p ~ 2M, and is 
continuous in the m —> limit, also at nonlinear level. 

The choice of Newtonian coordinates is allowed be- 
cause the coupling with matter preserves general covari- 
ance, and we also note that it enables us to show that the 
leading corrections start at 0(m 2 ) (unlike what would 
appear from the expansion of (|16p at small r, where ggg 
and g vv are not in the standard Newtonian form [loj). 

In the plain case (fT6|) of large-distance Yukawa falloff, 
a fairly strong bound on the scale to can be set from the 
existence of the largest gravitational bound states [2l| . 
which nowadays translates into the bound 1/m > 5 Mpc, 
or to < 10 -30 eV. This strong bound implies that the 
proposed models of massive gravity automatically pass 
the solar system tests, both at Newtonian and at PN 
level, given the above expressions (|T71) . 



As a phcnomcnological comment, we note that if Mi 
and M 2 are not of the same scale, e.g. if Mi 3> M 2 , and 
Ai < 0, then at intermediate distances the radial force is 
enhanced with the respect to the Newtonian case, leading 
to the flattening of rotation curves also in the absence of 
dark matter (before the Yukawa falloff at further larger 
distances). In this case, PN tests may become more rel- 
evant. We leave this scenario for a future investigation. 

Finally, let us briefly comment on cosmological so- 
lutions. In LI massive gravity, Friedmann-Robertson- 
Walker backgrounds simply do not exist [23l . |24| . Viable 
cosmologies do exist in the bigravity approach [24|, [25[ , 
where however cosmological perturbations are strongly 
coupled [2f|. In the present models it is likely that 
some extension will also be needed for cosmology (similar 
e.g. to |27f). Perturbations may instead turn out to be 
well-behaved, thanks to the absence of strong nonlineari- 
ties. We will address these issues in a separate work [l8j . 

Conclusions. We analyzed the general theories of mas- 
sive gravity with five propagating DoF and unbroken ro- 
tations, which can be candidates for healthy large dis- 
tance modifications of gravity. The conditions for these 
theories, formulated in [15j . were solved in general, in 
terms of three free functions. Then, as an example, we 
described a large class of models that can be analyzed 
explicitly, and have very interesting features. Namely: 
(a) Five DoF propagate both at linear and nonlinear 
level, (b) Gravitational waves are massive, including 
the three new graviton polarizations. These propagate 
with Lorentz-breaking dispersion relations, and may lead 
to striking signals in forthcoming gravitational wave ex- 
periments, (c) In static solutions, at short distances 
(tot « 1) the vDVZ discontinuity is absent, and there 
is no need to rely on nonlinear effects to pass the solar 
system tests of gravity; at large distances mr 3> 1 the 
gravitational potential is typically screened d la Yukawa, 
which is the desired behaviour of massive gravity, (d) In 
this case, one can set a limit on the graviton mass of 
m < 10 -30 eV. (e) As an effective theory the ultraviolet 
cutoff is A 2 = */ to Mp] , the highest possible in the ab- 
sence of a fundamental Higgs mechanism for gravity. Be- 
cause of the above bound, this corresponds to distances 
of ~ 10 mm, below which nonlinear and/or quan- 
tum corrections may lead to deviations from GR. This 
range is of interest for short-distance gravity probes [28| . 
(f ) Aside from this cutoff, the weak field approximation 
near a macroscopic source of mass M is valid down to the 
Schwarzschild radius Rg = 2M/Mp,, as in GR. Thus, a 
big advantage with respect to the Lorentz-invariant the- 
ories is that here modified gravity remains weak as GR 
in all the phenomcnologically interesting regimes. 

Of course to seriously compete with GR, massive grav- 
ity has to do more than just pass the solar system tests; 
it has to reproduce the many successes of GR, ranging 
from emission of gravitational waves from binary systems 
to CMB temperature fluctuations. Having a calculable 
weakly coupled theory at hand is of great help, if not a 
prerequisite, in carrying out the program. 
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Finally, we also remark that the conditions for 5 DoF, 
or uiq = at quadratic level, are not stable under quan- 
tum corrections, if one is willing to consider them. Vi- 
olations of these conditions spoil the counting of DoF, 



reintroducing the sixth mode, usually at the cutoff scale. 
This issue is shared among most massive gravity theo- 
ries, and calls for a proper symmetry protection for the 
framework described here. 
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